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ABSTRACT 



Finite-amplitude standing waves in air at ambient temperature 
contained in a rigid-walled rectangular cavity were experimentally 
investigated. The pressure waveform in the cavity was analyzed for 
harmonic content and compared to the theory of Coppens and Sanders 
as modified to include empirical losses. The waveforms for the 
nOO and OnO configurations were analyzed at strength parameters of 
about 0.5, 0.75, and 1.0, where the strength parameter is defined as 
MQi (y + l)/2. M is the Mach number of the fundamental, Q is the 
quality factor of the fundamental resonance, and y is the ratio of 
specific heats. It is found that the theoretical model accurately 
predicts the shape and approximates the amplitudes of the harmonic 
components of the finite-amplitude standing waves. 
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I. 



INTRODUCTION 



The purpose of this research is to make a detailed investigation 
of finite-amplitude standing waves in air at ambient temperatures in 
a rigid walled cavity, for comparison with the theory of Coppens and 
Sanders [l] modified to include observed energy losses. This in- 
vestigation is a continuation of the experimental works of Beech [2] 
and Winn [ 3D 3 with an extension to a rectangular cavity vice a long 
round tube. 

A rectangular cavity was chosen as it was a step toward the 
study of two and three dimensional waves, although this investi- 
gation concerns only one dimensional standing waves. In addition, 
it was thought that a rectangular geometry would produce over- 
tones which would be less harmonic than those for the more sym- 
metric geometry of the round tube. The dimensions of the cavity 
were chosen in an effort to space out the system resonances to 
avoid degeneracy or near -degeneracy . This was accomplished with 
some success, and thus allowed for the accurate location and anal- 
ysis of the first few overtones of each of the two modes studied. 
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II. 



BACKGROUND AND THEORY 



A finite amplitude wave is a solution of the non-linear wave 
equation that describes the real-life environment that confronts 
large amplitude acoustic waves. This finite-amplitude wave re- 
quires second order terms, or higher, to describe the waveform. 

It has been known for many years that hydrodynamic equations 
describing traveling waves in a non-diss ipat ive medium predict a 
change in the waveform. Kirchoff [4J and Lamb [ 5'J were among the 
early investigators in this area. In more recent years there have 
been many theories and methods of solution proposed to describe 
this distorted traveling waveform in a non-dissipative medium. 
Araoung these are the xvorks of Fay [6], and Blackstock [ 7 ] who pro- 
duced an extensive theoretical treatment of low amplitude non- 
linear waves, and made use of Fourier and power series 
appr oximat ions . 

When the problem was extended to dissipative mediums, it was 
shown very early [8] that waves in absorptive media exhibit finite- 
amplitude characteristics only at high strengths, different from 
waves in a non-dissipative media. More recently, an important con- 
tribution to the understanding of finite amplitude waves in a dis- 
sipative medium was presented by Fay L 9 3 • He, similar to a later 
work by Blackstock [lOj, formulated a Fourier series solution to 
describe the finite-amplitude waveform for a plane wave of infinite 
extent. More recent theoretical approaches to finite-amplitude 
effects have been of two basic types: an approximation solution, 
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such as that proposed by Fox and Wallace Cll], and a perturbation 
solution like that of Fay C 9 D • The approximation solution is based 
on assumed attenuation mechanisms and manner of overtone population. 
The perturbation solution makes use of perturbation methods and 
Fourier series representation of the waveform ClO, 12,13]. A per- 
turbation analysis which was extended to the sixth order by Keck 
and Beyer [13] is of particular importance to this investigation, 
as will be evident later . The solution obtained by Keck and Beyer 
is valid near the source, unlike previous perturbation solutions 
of finite-amplitude traveling waves in dissipative media. 

Although most investigations have been confined to traveling 
waves, there has been some investigations of standing waves. For 
example, Keller [l 4 l extracted solutions for standing waves in a 
closed tube. In that investigation, Keller assumed a non-dissi- 
pative medium which resulted in infinite amplitude waves at re- 
sonance, yet still obtained useful results for frequencies close 
to those of resonance. Coppens and Sanders Cl] developed an ex- 
tension of the Keck-Beyer perturbation approach, that include wall 
losses as predicted by the Rayleigh-Kir choff theory [l5l. This 
approach was carried to the tenth order by Ruff Cl6], at which 
point it was found that the mathematics become divergent. 

Coppens C 1 73 later expanded the theory and used a Fourier decom- 
position technique, which both Beech [ 3*] and Winn I2] experimentally 
investigated. It was found that the tubes that they investigated 
did not behave in the manner predicted by Ray leigh-Kir choff loss 
mechanisms. This investigation follows up the observation that 
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the tubes exhibited non-Rayleigh-Kir chof f behavior, by testing the 
theory using empirically determined losses. 

In the theory of Coppens and Sanders, a one dimensional, non- 
linear wave equation with dissipative terras which correspond to 
those encountered in the discription of plane standing waves in a 
rigid-walled cavity, is developed. The wave equation in this 
theory can be Fourier decomposed to take the form 
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and contains a dispersive term and a dissipative term, both based 
on the Rayleigh -Kir choff theory of wall losses C 1 5*1 9 and 

6 i = i/q i 

= particle displacement due to n^ harmonic 
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)^ + and CO^ are the half power frequencies above and below the 
fundamental . 

In the description of waves in cavities, Eq. 1 can be refor- 
mulated [l8] to introduce the actual cavity resonances and 
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generalized to include an arbitrary, empirical absorption coeffi- 
cient for each mode. The result is 
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= phase angle of n*^ harmonic component, where 0, = 0. 
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This equation can be manipulated into a general, time-independent 
form which, when evaluated at the rigid boundaries of the cavity, 
can be further simplified to produce the following form amenable to 
computer calculation: 



n = l 
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As can be seen from the above equations, it is important to know 
the Q(N) f s and OJ^’s of the system for the input of the theoretical 
predictions calculated by the computer. It is convenient to de- 
fine a quantity, E(n), that indicates the position of 00 relative 
to the classical harmonic frequencies, ntO^ . E(n) is defined by 

E ( n ) = ( w nr /( nw lr )) - 1 - 

The values of Q(n) and E(n) are directly determined from the 
infinitesimal amplitude analysis of the cavity and are the input 
quantities for the computer calculations. 
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III. EXPERIMENTAL CONSIDERATIONS 



A. APPARATUS 

A block diagram of the experimental apparatus is shown in 
Figure 1. Finite-amplitude waves were generated in air at ambient 
temperatures in a rectangular cavity with interior dimensions of 
6.97 cm. by 20.96 cm. by 30-28 cm. The walls of the cavity were 
constructed of O .985 in. thick aluminum plates, milled and screwed 
together to provide tight, right angle joints. These joints sealed 
by applying thin layers of silicon grease to adjoining surfaces 
before screwing them together. Acoustic waves were generated in 
the cavity by exciting the air with a plane-faced piston. The 
piston face gained access to the interior of the cavity by way of 
an orifice bored in one of the cavity walls. The piston was fitted 
with two O-rings and the piston orifice was bored such that the 
piston with 0-rings would make an airtight fit. It was found that 
the piston motion was restricted by the second 0-ring; therefore 
the piston was inserted into the orifice only as far as the first 
0-ring. Lubrication and further sealing of the pis ton -or if ice 
interface was provided by the same high vacuum silicon grease used 
to seal the wall joints. The piston was attached to an M-R Elec- 
tronics Model EA 1500 Exciter. This exciter was powered by two 
M-B Electronics Model 2120MB Power Amplifiers operating in parallel. 
The maximum attainable acceleration of the piston with this equip- 
ment was about 65 g T s in free air at 100 Hz. The frequency and 
power level of the piston were controlled by a General Radio Il6l-A 
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Coherent Decade Frequency Synthesizer. Frequencies can be selected 
to within less than 0.01 Hz with this synthesizer. The piston mo- 
tion was monitored by measuring the output of an Endeveco Model 2215 
Accelerometer mounted within the piston. Distortion in the piston 
motion, as determined by the accelerometer, set the lower limit to 
which reliable measurements of the non-linear distortion in the 
pressure waveform could be made. Measurement of the pressure wave- 
form in the cavity was made in the far corner from the piston 
source. Refer to Figure 2. The sound pressure levels and harmonic 
components of the pressure wave in the cavity were monitored by a 
1/4 inch diameter Bruel and Kjaer type 4136 Condenser Microphone. 
This microphone was mounted such that the active monitoring element 
of the microphone was flush with the interior of the cavity. The 
microphone was fitted with an 0-ring and seated in its orifice with 
the same high-vacuum silicon grease. The output signal of the 
microphone was amplified by a Bruel and Kjaer type 2801 Power 
Supply, This power supply with the microphone comprized the mi- 
crophone system. The microphone system was monitored by a Hewlett- 
Packard Model 400D Vacuum Tube Voltmeter. Another vacuum tube 
voltmeter of the same type was used to monitor the output of the 
accelerometer. A Hewlett-Packard Model 302A Wave Analyzer was 
used to measure harmonic distortion of both the accelerometer and 
microphone output signals. A Hewlett-Packard Model 521A Electronic 
Counter was used to monitor the filtered output of the wave ana- 
lyzer to assure that the proper frequency was being analyzed. 
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Figure 2. Rigid-Walled Cavity 



B. EQUIPMENT CALIBRATION 



Absolute calibration of the microphone system response was 
required to determine strength parameters. The microphone system 
included the microphone, its preamplifier, and the vacuum tube 
voltmeter connected to the output of the preamplifier. The micro- 
phone system response is defined by 

S = V /P 
m m 

where is the RMS voltage from the preamplifier as indicated on 
the VTVM and P is the RMS pressure amplitude at the microphone 
face. Calibration of the microphone system was accomplished with 
the use of a Breul and Kjaer Model 4220 Pistonphone which produces 
a 124 +_ 0.2 re volt/ubar acoustic signal at 250 Hz. The micro- 
phone system sensitivity was found to be 

S m = (1.46 + 0 .C 3 ) 10 “ 3 volt/ (Nt/m 2 ) , 
which is in close agreement with 

S = 1.50 + 0.03 x 10 _3 Volt/ ( Nt/m 2 ) 

m — v 

an< ^ -2 2 

S = 1.52 + 0.03 x 10 J Volt/ (Nt/m ) 

m — 

as determined by Beech [ 2 ] and Winn [ 3 ] respectively. The micro- 
phone system had a flat frequency response to within +_ 0.5 db 
from 80 Hz. to 20 kHz., as determined from the microphone response 
curve supplied by the manufacturer. With the microphone system 
calibrated, the wave analyzers were calibrated with respect to the 
microphone system VTVM. This was accomplished by removing the 
microphone and its preamplifier from the VTVM, and substituting 
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a voltage signal from the frequency synthesizer to both the VTVM 

and the wave analyzers simultaneously. 

The calibration of the accelerometer was accomplished by 

direct measurement of the piston face displacement at 100 Hz. and 

at various power settings with a stroboscopic light and a traveling 

microscope. The frequency of the strobe light was set very near 

100 Hz so that the beats of the piston/strobe light interaction 

could be observed. These beats, which occur at the difference 

frequency f , , - f . , , allowed the piston travel to be 

M * 1 strobe piston' 9 ^ 

easily measurable with the traveling microscope. Twelve deter- 
minations of the accelerometer sensitivity resulted in 

l p 

S = ( 2.36 0.04) * 10 Volt/(m/sec ). 

This compares with the value 

s = (2.35 + 0.03)- 10 -Zf Volt/(ra/sec 2 ) 
determined by Beech [33- 



C. MICROPHONE OUTPUT AND STRENGTH PARAMETER 

In order to investigate the pressure waveform in the cavity, 
as measured at the microphone, it was necessary to relate the 
strength parameter to observable quantities. The strength para- 
meter is defined as 



The peak pressure 



It is known from 



SP = MbQ . 



P^ is related 
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the microphone 



P = V/S 
rms ra 
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where V is the RMS voltage reading on the microphone system VTVM . 

The strength parameter can now be reformulated as 

SP = {2 VbQ 1 /(S m pc 2 ) , 

all of which are observable or calculable quantities. The calcu- 
lation or measurement of each of these quantities introduces error 
into the determination of the strength parameter; this error is 
estimated to be about 1.0 to 1.5 percent for both the modes 
investigated. 

D. FREQUENCY PARAMETER 

The frequency parameter is a quantity which indicates the 
position of the frequency being analysed relative to the fundamental 
resonance of the system. The frequency parameter is defined by 

FP = 2(W-W lr )Q 1 /C0 lr 

and can be shown to be related to the experimental frequency 
increment by 

FP = 2(C0-C0 lr )/(C0 1 + -C0 1 _) . 

E. SYSTEM ALIGNMENT 

The lower limit to which significant measurements could be 
taken was determined by the distortion introduced by the piston 
motion. In an effort to obtain meaningful data at the higher 
harmonics and for any harmonics at low strengths, great care was 
taken to reduce the distortion introduced into the cavity. The 
distortion of the piston motion when free of the cavity was below 
the general electronic noise level of the accelerometer system and 
thus was insignificant. The bulk of the distortion came from 
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misalignment of the piston and cavity orifice. The procedure used 
to obtain maximum alignment of the piston and thus minimize the 
distortion was as follows: The output signal of the accelerometer 
system, which measured the piston motion, was displayed on the 
dual trace oscilloscope and this waveform was compared to a pure 
sine wave. The piston and exciter apparatus, which was mounted in 
an adjustable framework, was positioned relative to the cavity 
orifice until the accelerometer waveform was as clean and undis- 
torted as possible. In order to minimize any effects which may 
have arisen from vibration coupling of the cavity and exciter or 
cavity and table, the cavity was mechanically isolated from both 
the exciter and the table. This was accomplished by mounting both 
the cavity and the piston-exciter apparatus on platforms con- 
structed of alternate slabs of wood and rubber. This restricted 
the interaction of the cavity with the piston to the piston ori- 
fice only. With proper alignment, it was possible to reproduce 
conditions where the distortion in the input waveform to the 
cavity was less than 0.2% second harmonic and 0.3% for third har- 
monic, as measured by the H-P 302 a wave analyzer. Distortion 
values for the higher harmonics were below the general ambient 
noise of the total apparatus, and thus were ignored. 
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IV. PROCEDURES 



A. PRERUN PROCEDURES 

Before any meaningful data could be taken, it was important to 
determine the amount of vibration transferred through the cavity 
walls and picked up by the microphone. To make this determination, 
the top of the cavity was removed, thus relieving the pressure seal, 
and a wood plate added between the piston and microphone to isolate 
the microphone from the pressure waves induced in the air by the 
piston. As described in the System Alignment section, the cavity 
and the exciter -pis ton assemblies were isolated by a floating 
platform technique. The vibration picked up by the microphone was 
measured at frequencies of interest and produced microphone system 

-5 -4 

responses of 8*10 to 2*10 volts. In all cases these values 
amounted to distortion well below 0.3%, which is the maximum 
amount of piston motion distortion. 

It was noted by Beech [ 3 ] that temperature changes produced 
significant changes in the resonant frequency of the tube he in- 
vestigated. The temperature effect on the cavity investigated 
here was very pronounced. This presented many problems as the 
frequency of fundamental resonance was the reference for all de- 
terminations of frequency parameter. Therefore it was necessary 
to know the value of the fundamental resonance at all times. To 
minimize this shift in the resonance frequency, the entire ap- 
paratus was allowed to warm up for a period of one to four hours, 
depending on the strength parameter to be investigated, prior to 
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taking data. This procedure did not totally relieve the effect of 
the frequency shift so another calibration method was employed. 

Many determinations of the fundamental frequency were made during 
the taking of data, thus producing a calibration curve, see 
Figure 3? by which the frequency parameter values could be 
cor r ected . 

Immediately prior to the taking of each data set, the align- 
ment of the piston and cavity was checked as outlined in the system 
alignment section of Experimental Considerations. From this point 
onward, the procedure for data collection was dictated by the 
quantities being measured. There was one experimental method to 
measure the Q T s and E f s of the system in the linear region, and 
another to measure the harmonic content of the non-linear finite 
amplitude waves. 

B. MEASUREMENT OF Q(N) AND E( N ) IN THE LINEAR REGION 

Data collection in this region required very little power, 
but required very fine control of the power setting. This control 
was obtained by inserting a voltage dividing circuit between the 
frequency synthesizer and the power amplifiers. This voltage di- 
vider reduced the available power by a factor of ten and increased 
the sensitivity of the power control by the same factor of ten. 

With this arrangement, a satisfactory fine tuning power control 
was obtained at low strengths. 

With the power set at some nominal low value, the frequency 
region about the resonant or harmonic pressure peak of interest 
was scanned until a maximum in pressure was indicated by the 
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Figure 3. Fundamental Frequency Calibration Curve 



microphone system response. The input power, determined in the 
linear, low strength, region, and kept constant during the fol- 
lowing measurements. The frequency of the driving piston was 
changed above and below the frequency of maximum pressure, until 
the pressure was -3 dB of the maximum pressure value. This deter- 
mination of the -3 dB points was made with the 3°2A wave analyzer , 
which measured the pressure amplitude of the cavity by way of the 
microphone system. 

It was found that the first seven overtones of the nOO mode 
and the first five overtones for the OnO mode were undisturbed 
by overlapping Q curves. The value of Q for each response curve 
was determined from the frequency difference of the -3 dB points, 
divided by the center point frequency between the two -3 dB 
points. This center point was assumed to be a good estimate of 
the peak pressure frequency, CO^ , for the response curve in 
question. It was with these values for the system resonance fre- 
quencies that the E(n) T s were determined. 

C. MEASUREMENT OF THE HARMONIC CONTENT OF FINITE AMPLITUDE 

STANDING WAVES 

In order to obtain finite amplitude waves in the cavity, it 
was necessary to go to much higher strengths than were required 
for the measurement of the Q f s and E f s. With the driving fre- 
quency set to some value near the fundamental, and the peak pres- 
sure component of the fundamental held constant, the strength of 
each harmonic was measured by one of the wave analyzers. The 
other wave analyzer monitored the fundamental pressure , which 
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was kept constant by varying the power setting of the synthesizer 
as necessary. Values of the frequency parameter were chosen such 
that the distortion curves, Figures 4 through 9 , could be deter- 
mined in the areas of maximum distortion. The order in which 
measurements were taken for each distortion curve was randomized 
to avoid systematic error. Strength parameters of about 0 . 5 ? 0 . 75 , 
and 1.0 were investigated for each mode. Strengths below 0.5 did 
not produce high enough distortion to be measured with the equip- 
ment available. Strengths above 1.0 were not accessable due to 
limitations in power. 
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V. RESULTS 



Tables I and II are the data for the linear, infinitesimal 
amplitude region. The information on the empirical losses is con- 
tained in the Q ! s and E T s. These are the values used in the com- 
puter program to predict the harmonic distortion curves in Figures 
4 through 9. The theoretical curves in these figures are plotted 
along with the experimentally measured values of harmonic dis- 
tortion. The theoretical predicted values were generated for fre- 
quency parameter intervals of 0.2 and were plotted as solid curves 
The experimentally measured values were plotted as dots where they 
occured. Data were taken, and theoretical predictions made, for 
strength parameters of 0.513, 0.782, and 1.00 for the nOO mode and 
0.597, 0.780, and 1.00 for the OnO mode. 

The point of intersection of the error bars on the P^/P^ theo “ 

retical curve in Figure 4 corresponds to the value of EP = 2Q E(2) 

CO 

2r 

where CO = — - — . It is important to note that the n = 2 distortion 
curves peak when the system is driven at this frequency. That is, 
when the driving frequency CO is equal to the classical fundamental 
of.CO^r* there is maximum content of P Q . The point where this oc- 
curs for each P^/P^ curve indicated by arrows in Figures 5 
through 9 and by crossed error bars in Figure 4 . The error in the 
determination of these points comes mainly from the error in E(2) 
that appears in the above equation for FP . The error in is of 
much less significance as is indicated in Table III. The error in 
the P^/P^curves in Figures 5 through 9 are approximately the same 
as those indicated in Figure 4 * 
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Table 1 



nOO Mode 



n 


f.(Hz) 


f + (Hz) 


Af (Hz) 


f (Hz.) 


Q 


E x 10 3 


i 


677.80 


679.70 


1.90 


678.75 


357.24 


0 


2 


1365.74 


1368.49 


2.75 


1367.12 


497.13 


7.0 


3 


2059.17 


2062.84 


3.67 


2061.01 


561.58 


12.0 


4 


2787.37 


2792.73 


5.36 


2790.05 


520.53 


27.5 


5 


3423.00 


3427.80 


4.80 


3425.40 


713.63 


9.1 


6 


4126.41 


4131.91 


5.50 


4129.16 


750.76 


13.6 


7 


4803.92 


4809.93 


6.01 


4806.93 


799.82 


11.4 


1 


678.02 


679.95 


1.93 


678.99 


351.81 


0 



OnO 1 
n 


Mode 

f_(Hz) 


f + (Hz) 


Table II 
Af(Hz) f n (Hz) 


Q 


E x 10 3 


1 


813.99 


816.08 


2 .09 


815. c4 


389.97 


0 


2 


1633.16 


1636.39 


3.23 


1634.78 


506 . 12 


2.8 


3 


2415.11 


2420.57 


5.46 


2417.84 


442.83 


-11.3 


4 


3143.74 


3150.10 


6.36 


3146.92 


494.80 


-34.9 


5 


39H.05 


3920.35 


9.30 


3915.70 


423.62 


-39.4 


1 


814.23 


816.34 


2.11 


815.29 


386.39 


0 
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Frequency Parameter (FP) 
nOO Mode SP = 0.513 

Figure 4 
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Frequency Parameter (FP) 
nOO Mode SP = 0.782 

Figure 5 
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Figure 7 
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Table III contains the results of five separate determinations 
of the Q T s and E f s taken in the infinitesimal amplitude region. 
These measurements were made to provide an estimate of the error 
in the theoretical curves. In all cases the values used to gene- 
rate the curves were very close to the average values determined 
later and the regeneration of the curves was therefore unwarranted. 
Table III contains both sets of values for direct comparison. 
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Table III 



nOO Mode 



Average Values 



Values Used to 
Determine Curves 



n Q(n) 



E(n) ’ 10 3 



Q(n) 



E(n) • 10 3 



1 


358.2 +0.4 




357.2 




2 


498.5 + 4.7 


6.8 +0.1 


497.1 


7.0 


3 


559.8 +0.8 


11.6 + 0.2 


561.6 


12.0 



OnO Mode 



1 


359.8 


+ 5-3 




390.0 




2 


515.7 


+ 5.1 


2.6 + 0.1 


509.1 


2.8 


3 


459.2 


+ 7.5 


-11.3 + 0.1 


442 .8 


-11.3 
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VI. CONCLUSIONS 



It appears that the present theoretical model successfully 
predicts the major features of the harmonic content for finite- 
amplitude standing waves in a rigid-walled cavity. In all cases 
investigated, the shape and general magnitude were predicted; see 
Figures 4 through 9. 

It is interesting to note the position of each theoretical 
curve relative to its corresponding experimental curve. For the 
nOO mode, the theoretically predicted values are consistently 
above the measured values. The situation for the OnO mode is the 
same for P /P^ but opposite for P^/P^. A possible source of this 
difference is the large difference in the geometry encountered by 
each mode. The nOO mode is a standing wave whose wave vector is 
perpendicular to the wall with the microphone orifice, which must 
be considered a very small perturbation due to its small size. On 
the other hand, the OnO mode is a standing wave with its wave 
vector perpendicular to the piston orifice, which is a large per- 
turbation in the geometry. In addition, the fact that the wave 
generating surface, the piston face, is normal to the nOO wave 
vector and perpendicular to the OnO vector could cause a con- 
sistant difference in the behavior of the modes. It is certainly 
unclear to this investigator just what the specifics are that 
explain the trends in the theory/experiment relationship. It is 
clear however, that the theory of Coppens and Sanders as modified 
to include empirical losses, does predict the essence of the 
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behavior of finite-amplitude standing waves in the rigid-walled 
cavity investigated. 
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APPENDIX A 



In the investigation by Winn [ 2 ], it was found that the tubes 
analyzed did not behave in a manner consistent with Rayleigh- 
Kirchoff wall losses. In an effort to determine if the interior 
finish of the tubes was a critical factor that determined the loss 
mechanisms, the six foot tube was reanalyzed after refinishing the 
interior to a highly polished state. These results were compared 
to the results for the tube in a smooth but not polished condition 
as determined by Winn. 

It was found that the absorption coefficients, defined as 

CO - CO 

OL = — — — 

n 2C 

o 

were not significantly changed from the values found by Winn, and 
both sets of values obeyed the Ray leigh-Kir chof f predictions 
indicated by the solid line in Figure 10. The behavior of the 
overtones (for the polished tube) was found to be significantly 
different (from the overtone behavior in the unpolished tube), 
and was significantly different from the Ray leigh -K ir chof f pre- 
dictions as indicated by a solid line in Figure 11. 
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Figure 10. Comparison of Attenuation Coefficients 
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Figure 11. Comparison of E(n)s 
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